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Abstract 

We investigate the interaction between spherical magnetic nanoparticles which present either a single 
domain or a vortex structure. First the magnetic structure of a uniaxial soft sphere is revisited, and then the 
interaction energy is calculated from a micromagnetic simulation. In the vortex regime the orientation of 
the vortex relative to the easy axis depends on both the particle size and the anisotropy constant. We show 
that the leading term of the interaction is the dipolar interaction energy between the magnetic moments. 
For particles presenting a vortex structure, we show that the polarization due to the dipolar field must be 
included. The parameters entering in the dipolar interaction are deduced from the magnetic behavior of the 
isolated particle. 
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1 Introduction 

The magnetic behavior of nanometric particles either isolated or in nanostructured bulk materials is now quite 
well undertood both from experiments or numerical calculations [1] [5] . With the growing diversity of systems 
made of such nano-objects as building blocks either as 2D or 3D systems in non magnetic environment or in 
colloidal suspensions as ferrofluids, a precise knowledge of interparticles interaction is needed. This has already 
been done at least partially in a variety of systems, such as nanograins OH] nanorings [S] or flat nanodots [S1I3- 
with a predominant attention paid on short range effects. A lot of work remains to be done in this field and 
especially for spherical particles. In particular it seems important to developp models for the long ranged and 
anisotropic dipolar interaction. A further interpretation of experimental results as those of [5] necessitates such 
models. 

2 Magnetization structure and hysteresis 

We consider a spherical particle characterized by a radius R ranging from lOnm to Abnm, an exchange constant, 
Ax = 1.10^^^ J/m and a saturation magnetization Jg = IT corresponding roughly to Permalloy. The anisotropy 
constant of the uniaxial magneto-crystalline energy is varied between Ki = and 7.10^J/m^. The numerical 
calculations are performed with the the framework of micromagnetism with the code MAGPAR [9] based on a 
FEM method. The particle volume is Vs and hatted letters denote unit vectors. Small particles, up to roughly 
R = 20nm for our parameters, are uniformly magnetized as single magnetic domains with square hysteresis 
curve. When the particle radius increases beyond roughly R = 20nm, a vortex structure is obtained, and the 
local magnetization profile, TO(r), is decomposed in its cylindrical components using the vortex axis, say i), as the 
cylindrical axis. At zero or small values of the external field, the direction of the vortex relative to the easy axis, 
a, depends on both the value of R and of iCi, a behavior already obtained in the case of the cubic anisotropy 
[TU] . We find, as expected, w ± a either for R > Rth or Ki > Ki^th at constant Ki or R respectively, Rth and 
Ki^th being some threshold values. [Rth = 26nm for Ki = 3. 10^ J/m^ and Ki^th = 2.10^^J/m'^ for R = 45nm). 
The orientation of i) relative to a affects strongly the M{Hex) curve. Whith an external field along a, w _L a leads 
to a magnetization curve qualitatively similar to that of a flat nanodot with an in plane external field [2] : 
threre is no remanence and M{Hex) presents two lobes. For spherical particles these lobes are associtated to 
the rotation of the vortex core in the direction of the field prior to its annihilation. The magnetization normal 
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Figure 1: Magnetization parallel and normal to the field, hex — a. (solid line) or v (dashed line). R = Abnm, 
Ki = i.W^J/m. 

to the field takes a non zero constant value corresponding to the vortex core magnetization, when the variation 
of M \\ results from a shift of the vortex normal to the field, this value coincides with the remanence obtained 
with hex ~ i>- This behavior is displayed on figure [TJ 

As is generally obtained in nanodots or spherical particles [TT], O [I3| the magnetization M in the direction 
of the external field depends nearly linearly on the field in the vicinity of Hex = and away from swhitching 
points where the vortex reverses as a whole. Such a linear behavior is observed both when hex = v or hex 
_L V. This means that the susceptibility x defined as = X does not depend on the field to a very 

good approximation. This can be exploited to obtain the variation of the total energy with respect to the 
external field. We analyse the variation of the magnetization, AM as the polarization of the sphere induced 
by the field. Starting from AM{Hex) — X^ex, we deduce AAI{Hex) from the energy, E{Hex) by writting 
an equilibrium equation '^^q^^P = 0. The total energy depends explicitly on Hex through the Zeeman term, 
Ez = — ^{)Hex{'m[Q)v .hex + AM), where we have expressed the permanent magnetization in the absence of the 
field as M{Hex = 0) = m(O)i) and m(0) denotes the magnitude of the vortex core magnetization in the absence 
of the field. Then we get 
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Therefore the variation of the total energy is 

E{Hex)~E{0) 



{Edm + Ex + Ea) — floHe 
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AM' + Ez{Hex) 
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where we have used AAf (Hex) ~ xHex- Mo is the vacuum permeability and Ez the Zeeman energy. Notice that 
([T|) is exact while ^ holds only in the case of a linear dependence of AM {Hex) with respect to Hex- The first 
term of the r.h.s. of ^ is the polarization energy of the sphere [T3] and corresponds to the energy cost for the 
reorientation of the magnetization inside the sphere. It can be written equivalently as (fj,QAAIHex)/2 

3 Interaction between magnetic nanospheres 

The interaction energy between two magnetic nanoparticles in terms of the interparticle distance, ri2, is defined 
from the total energy of the two spheres system 
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where (1,2) is a short notation for the orientation and location variables of the particles. We expect to get a 
form dictated by the dipolar interaction between the magnetic moments of the approaching spheres wich reads 

£^d,p = ^^£^dii2(mi,m2,fi2) (4) 

(iii2(mi,m2,^i2) = mi.7fi2 - 3(mi.fi2)(TO2-»^i2) (5) 

where mi are the magnitude of the magnetic moments. For single domain particles, = MsV^ where is 
the saturation magnetization of the particles, and the orientations rhi result from the minimum of dii2 given 
in ([5]). For particles without magnetocristalline anisotropy, this gives: rhi = m2 = fi2 and dii2 = -2. For non 
zero magnetocristalline energy on both particles with easy axes 0,^, the orientations rhi result of the interplay 
between the anisotropy energy tending to align rhi on and the energy ^ tending to minimize the angular 
function If Ki takes a non vanishing value only in one particle say i = I and is large enough to impose rhi = 
0,1, 7712 must orient in the dipolar field due to particle 1 i.e. in such a way that dii2 = m2-ai — 3(TO2.fi2)(ai.fi2) 
is minimum. The whole behavior outlined above is very well reproduced by the full micromagnetic calculation 
which demonstrates the dipolar nature of the interaction between single domain particles. 

In the case of particles large enough to present a vortex structure, the orientations of the effective moments 
of the particles arc the vortex directions, iii, and the values of the moments correspond to the vortex cores 
magnetizations. We introduce the coefficient = mi/{MsVs). The value taken by a depends on both the 
characteristics of the particle and ri2 and c?ii2 through the polarization of the particle by the dipolar field of 
the second one. (q;(1,2) = a(ri2,(iii2) and ao ~ a{oo)). A simple approximation for the interaction energy 
is built by considering that each particle is in the dipolar field of the other. We have to take into account 
two contributions. The first one is given by and the second one is twice the polarization energy of one 
sphere in the field of the second one. The second contribution has been introduced in ^ for one particle in 
a constant external field. The role of m(o) is played by MsVsCto while the induced moment in the direction of 
the dipolar field is p = phdip = X H dip {r 12) h dip. We consider the case where the vortex Vi is free to orient in 
the direction of the dipolar field due to particle j ^ i. This corresponds to either the absence of anisotropy 
or particles large enough for the vortex to be normal to the easy axis and oi = 02. In this case we have 
p ~ {a{l, 2) — ao))MsVshdip = Aq;(1, 2)MsVsV and adding twice the first term of (U) to the dipolar energy we 
get 

S,„t(l,2) = MM^c,o(aj, + Aa(l,2))dn2 (6) 

which coincides with the interaction energy between polar polarizable hard spheres. Then we have to calculate 
Aa(l,2); for the simple case of two particles, introducing u = x/{4:-7tR'^) we get 

Aa{l,2) = a,—-g^P^^-— (7) 
(ri2/i?)'^ + udii2 

A typical example of the interaction in terms of the distance ri2 for dii2 ~ —2 is displayed on figure [H 
Similar results are obtained for other particles. The full calculation is compared to the analytical model given 
by ([6|). We consider two levels of approximation: either the simple dipolar interaction, where the polarization 
energy is neglected {u = 0) and the interaction corresponding to the dipolar polarizable spheres. Except in the 
vicinity of contact, (ri2 — 2R) equ. (O reproduces quite well Ei„tii,2). The parameters needed, ao and x a-re 
deduced from the magnetization curve of the isolated particle. We can therefore conclude that the interaction 
between spherical soft magnetic particles is of dipolar nature, and that the polarizability must be included when 
the particles present a vortex structure. 



4 Concluding remarks 

We can conclude that the leading interaction between soft magnetic nanoparticles corresponds as expected to 
the magnetic dipolar interaction. Then two situations must be distinguished according to the particle size: the 
interaction corresponds to the total magnetic moment without any polarizability contribution in the case of small 
particles which are in a single domain state when isolated, while a polarizability term must be included when 
the particles present a vortex structure. In this latter case, the value of the permanent moment, represented 
by the coefficient ao in the present work, is no more trivial, but can be determined as well as the relevant 
susceptibility x from the isolated particle properties. We emphasize that here, the solvation od the equation for 
the variation of Aa in terms of ri2 and ^112(1, 2) is quite obvious since there is only two particles, but this should 
be not the case for an assembly including a large number of particles where the interaction will present a n-body 
character. Finaly, because of the the dipolar nature of the leading interaction, we can confirm the behavior 
observed experimentaly in [S] since the vortex of neighboring particles are expected to allign themselves. 
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Figure 2: Interaction energy per unit volume between two spheres. R ~ 35nm; Ki ~ 0; dii2 = —2. Open 
circles: full calculation, from Q (the thin line is a guide to the eye); solid line: equ.®; dashed line: simple 
dipolar approximation, u ~ 0. 
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